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ABSTRACT 

A  reliable  system  for  performing  a  Fourier  analysis  of  a 
reactor  noise  signal  has  been  constructed  and  utilized.   Both 
auto-  and  cross-correlation  techniques  were  used  to  study 
noise  signals  from  the  MITR  hohlraum  and  the  exponential 
lattice.   As  predicted  by  theory*  the  efficiency  of  the  detec- 
tors used  was  insufficient  to  provide  positive  results.  The 
lattice  transfer  function  was  electronically  simulated  and 
analyzed  by  the  system  confirming  the  reliability  of  the 
equipment . 

Several  methods  are  suggested  for  increasing  detector 
efficiencies  to  provide  signal  levels  which  will  permit 
extraction  of  information  about  the  kinetic  behavior  of  the 
test  lattice* 
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CHAPTER  1 
INTRODUCTION 

1 ol  Background 

In  studying  the  neutron  characteristics  of  any  nuclear 
reactor  core  assembly  It  is  desirable  that  the  research  be 
done  as  economically  as  possible ^  and  preferably  that 
several  independent  methods  be  used  to  determine  the 
desired  information.  The  MIT  Lattice  Research  Project  has 
for  some  time  studied  various  reactor  parameters  in  a  heavy 
water-moderated  exponential  facility  located  adjacent  to 
the  MIT  Reactor «v  -"     See  Figure  1,1  for  a  cut-away  view 
of  the  MIT  Lattice  Research  Facility „ 

In  recent  years ,  considerable  use  has  been  made  of 
the  pulsed  neutron  technique  in  determining  k„  for  test 
lattices  in  order  to  provide  data  supplementing  the  more 
common  static  parameter  measurements  <S    '     It  is  the  purpose 
of  this  thesis  to  explore  the  utility  of  the  noise  analysis 
technique  as  a  means  of  complementing  the  pulsed  neutron 
and  other  methods  of  measuring  the  various  parameters 
measured  in  the  MIT  heavy  water  lattices  „ 

1»2  Noise 


The  term  "noise M  refers  to  the  random  or  stochastic 
fluctuations  in  the  output  signal  from  a  system 0  Reactor 
noise  is  thus  the  random  fluctuation  of  power  level  due 
to  the  discrete  nature  of  the  multiplication  process  „ 
These  fluctuations  are  generally  small ,  and  usually  only 
about  one  percent  of  the  total  power  «  In  order  to 
observe  this  noise,  a  detector  must  be  used,  and  often 
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Fig.  1.1  Cut-away  View  of  MIT  Lattice  Research  Facility 
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the  terra  reactor  noise  is  used  to  describe  the  fluctuations 

of  the  output  signal  of  the  detector •  It  is  important  to 
remember,  however,  that  the  detector  only  seej/s  a  token 
sample  of  these  fluctuations  aad  that,  owing  to  the  finite 
size  of  the  detector,  the  spatial-dependence  of  the  noise 

measurement  is  blurred „  Only  recently  ha3  the  space-depen- 

(7) 
dence  of  reactor  noise  been  closely  examined »*  '  Throughout 

this  paper,  however,  space-independence  will  be  assumed 

since  this  will  prove  adequate  fo  b* 

1o3  Summary  of  Theoretical  Considerations 

Since  a  reactor  depends  upon  neutron  multiplication  for 
its  operation,  it  is  reasonable  to  expect  that  among  the 
various  neutrons  which  are  registered  by  a  detector,  many  which 
have  a  common  ancestry  will  produce  pulses  superimposed  upon 
the  doCo  level  of  the  signal .  The  number  of  such  related 
pulses  will  depend  upon  the  detector  efficiency  as  will  be 
seen  later .  While  the  detector  d0c«  level  is  itself  main- 
tained by  the  accumulation  of  many  pulses,  we  are  unable  to 
extract  any  information  involving  the  discreteness  of  the 
multiplication  process  from  it,  and  it  thus  holds  no  interest 
for  us  here  8  Thus  one  can  instead  focus  attention  upon  a 
fluctuation  in  the  signal  due  to  the  daughter  neutron  of  a 
given  fission,  followed  a  certain  At  later  by  another 
impulse  due  to  a  grandaughter  neutron  of  the  same  fission. 
Since  this  time  interval  At  is  a  func     of  the  neutron 
lifetime  in  a  given  system,  it  should  be  poss     to  extract 
information  about  the  multiplying  ]         of  the  medium 


by  submitting  the  aoCo  (or  noise)  portion  of  the  detector 
signal  to  analysis  „  One  method  of  analysis  is  to  determine 
the  autocorrelation  function*  which  is  defined  as  the  expec- 
tation value  of  the  product  of  the  signal  at  a  time  tj,  and 
the  signal  at  a  time  t^s 

0nn(t)   =E[n«tI)n(t,ij   ,  (1.1) 

where  for  simplici^      signal  is  taken  to  be  the  neutron 
density  n(t),  and  T  =  I  evelopment  of 

Bierman,  et  al.,'8^  one  can  evaluate  0nn(Z)   by  first  con- 
sidering the  standard  reactor  kinetics  equations  in  their 
space -independent,  linear ,  time -invariant  forms 

dn  =  kU^M  n  +  s  x.C,  +  S  (1.2) 

clC,  p.  . 

dt      i  i    £ 

Since  only  prompt  neutron  behavior  is  of  importance* 
these  equations  may  be  put  in  the  forms 

^+  an  =k|=n0  +  S  -  f(t)  (1.4) 

where  .  „AdjfcftL  ,  the  prompt  neutron  decay  constant,  and 
the  inhomogeneous  terms  are  considered  as  the  forcing  func- 
tion f(t). 

The  solution  to  Eq.  (1.4)  may  be  written  as  a  convolu- 
tion integrals 

/CO  if00 

n(t)  =     W(t-t)f(t)dt  :    W(t)f{t-r)    ,(1.5) 
J  -co  ) 


where  w"(t~t)  is  the  impulse  response  function  of  the  system. 
This  function  is  obtained  by  normalizing  the  homogeneous 
solution  to  Eq0  (1.4)  and  is  found  to  beg 

W(t-£)  =  I.  (1.6) 

o  t  <  tr 

One  may  evaluate  Eq.  ("*  1)     igs 

fijv  -I e'at-  (1°7) 

If  one  were  to  assume  that  the  forcing  function  were 
instead  a  burst  of  neutrons ,  as  in  a  pulsed  neutron  experi- 
ment, f(t)  could  be  approximated  by  a  Dirac  delta  6(t)  giving 
as  a  solution  to  Eqa  (1.4)s 

n(t)  =  ae~at  o 

At  this  point,  the  complimentary  nature  of  the  pulsed 
neutron  and  noise  analysis  experiments  becomes  evident,  for 
each  provides  an  independent  method  for  determining  the 
prompt  neutron  decay  constant  a. 

Various  noise  analysis  methods  are  available  for  deter= 
mining  a  by  studying  the  a  0c  .  detector  signal  in  the  time 
domain,^)  however,  by  transforming  Eqa  (1.5)  into  the  fre- 
quency domain  by  the  use  of  Fourier  transforms,  a  method  of 
analysis  which  involves  less  data  processing  can  be  employed* 

The  Fourier  transform  of  the  autocorrelation  function 
is  called  the  power-spectral=density  function  and  is  given 
bys 

J   -00 
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which  when  combined  with  EqQ  (17), 

">J-*fo]  (1-9) 

The  spectral  shape  of  Eq.  (1»9)  is  that  of  a  low  pass 
filter  with  a  breakpoint  at  a     .  Thus,  if  a  Fourier  analy- 
sis of  the  detector  signal  is  performed,  it  should  be  possi- 
ble to  determine  a,  the  p  ;ay  constant „ 

lc4  Detector  Efficienc; 

In  order  to  co  neutron 

and  noise  measurers  '  noted  that  both 

experiments  are  hampered  by       spurious  signal  levels  •  In 
the  pulsed  neutron  experiment     3e  take  the  form  of  residual 
neutron  levels  due  to  the  neutron  sources  in  the  system,  and 
can  be  overcome  by  increasing  pulse  intensity „  In  noise 
experiments,  as  we  shall  see  in  this  section,  the  detection 
process  Itself  provides  a  white  background  power  spectrum, 
which  can  only  be  overcome  by  increasing  detector  efficiency 
or  by  performing  cross-correlation  expei    fes0  These  limita- 
tions can  be  described  by  including  the  background  as  an 
additive  constant  to  both  experiments s 

N(t)  =Ap  +  Bpe-ot  (1.10) 

(loll) 


1  +-J2/o2 


where  now  a  figure  of  merit  for  each  system  can  be  defined  ass 


B 

A  +  B 
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The  utility  and  feasibility  of  noise  experiments  thus  depends 

upon  the  magnitude  of  the  figure  nt   mp-nit. 

By  oursuing  the  analvtical  method  of  Gohn*  '  ' .  and  extend- 

ing  it  to  the  case  of  a  subcritical  reactor,  the  effects  of 

detector  efficiency  can  be  studied „  The  basis  of  this  method 

rests  on  what  is  known  as  the  Schottky  formula,  which  was 

originally  developed  to  describe  the  noise  in  a  temperature 

limited  electronic  diode 0  It  iss 

<|XI2>  =  2e2m  ,  (1.12) 

o 
where  <fl|  >  is  the  spectral  density  of  the  diode  current 

p 

noise  in  amps  sec,  e  is  the  charge  carried  by  each  electron 

in  coulombs,  and  m  is  the  average  number  of  electrons  flowing 
per  second  a  Cohn  adapts  this  formula  to  the  reactor  case 
giving § 

(!S0l2)=2Zq^ilJ  (lol3) 

where  OS  l2>  is  the  spectral  density  of  the  noise  equiva- 
o 

lent  source  in  neutrons2  sec  ,  q^  is  the  net  number  of  neu- 
trons produced  in  a  nuclear  process  of  type  i,  and  m±   is  the 
average  number  of  nuclear  processes  of  type  i  occurring  per 
second  in  the  reactor .  As  show  in  Appendix  A  this  formula- 
tion results  in  the  following  expression  for  the  noise  equiva- 
lent source  in  a  subcritical  assemblys 


OS0!2>  =  s 


.  k  .  iL+iL^-  2  v  +  1)    (1-14) 

This  source  is  transformed  by     reactor  scarce  transfer 
function  T(«u>)  into  the  neutron  s:  sitys 


<|n(ej)l2>  =  \T(co)|2<|S  (lol5) 

By  considering  the  power  spectral  density  of  the  detec- 
tor current  due  to  both  the  correlated  and  uncorrelated  neu- 
trons, it  can  be  shown  that  the  total  urrent  auto 
power-spectral-density  1st 


<UaM'2>    =2Q2€P0jl  +  r~f 


°kD+  J£{7^2V+1) 


'  3 — g  > 


(1.16) 


where 
€ 

F 


o 


k 


P 


V 

Q 


«*.» 


_  fissions  detected 

m  detector  efficiency  -  fi8Si0ns  occurring 

k  n 
=  fission  rate  =  ■— p 

=  k(i-P) 

=  average  number  of  neutrons  emitted  per  fission 
=  charge  collected  per  fission  detected 
1 


Thus,  for  an  autocorrelation  measurement,  the  detector 
efficiency  must  be  such  that  the  right-hand  term  is  greater 
than  unity,  or  for  <^«  as 


(1~V 


>      V  Br,-  ■■— — 

L     v     p  J 


(1.17) 


Using  the  accepted  values(9J  of  2*5  for  the  factor  y  and 
7o4?  for     $   ,   this  becomes s 

k^50-~V  _  Cia8) 

kp[|  -  3o5  kp  +  3o47  kp*] 


el 


Figure  lo2  displays  this  ■  Oo90o 

Since  detector  efficiencies  above  lCT^  are  very  difficult 
to  obtain  without  significantly  disturbing  the  system  under 
study 9   a  major  sturabll-  he  use  of  noise  analysis 

in  far-subcritical  systems  {       becomes  obvious „  Some  of 
this  problem  could  be  alleviated  if  the  uncorrected  noise  is 
filtered  out „  By  us     wo  detectors  and  multiplying  the  a  0c . 
components  of  the!  ose  pulses  which  did  not  coin- 

cide would  tend  to  Cancel  out  where  h  did  coincide 

would  reinforce  each  other  and.,  over  a  period  of  time*  the 
correlated  portion  of  the  spectrum  would  tend  to  foe  emphasized „ 
This  method  of  cross-correlation  was  used  by  Seifritz, 
Stegemann,  and  Vath  in  making  measurements  on  the  coupled  fast- 
thermal  reactor  STARK o^11^  Their  results  indicated  that  with 
their  equipment*  a  ratio  of  correlated  to  uncorrected  noise 
of  Od  was  sufficient  to  give  meaningful  results  if  cross- 
correlation  techniques  were  used,         a  factor  of  twenty 
improvement  over  the  rati     2  which  they  needed  in  auto- 
correlation experiments.  Of  course,  the  relative  advantage  of 
cross-correlation  over  auto-correlation  experiments  depends 
upon  the  individual  analysis  i    uents  and  how  the  signal  is 
analyzed .  Figure  1  1  shows  the     -correlated  detector 
efficiencies  required  based  on  their  experience ,  Notice  that 
detector  efficiencies  of  close  to  one  percent  are  still 
required  for  a  k  of  0»T  to  Qu8,  which  is  the  range  of  pri- 
mary interest  in  MIT  lattice  c. 

1°5  Purpose  of  Present  Work 

This  thesis  is  concerned,  the,  assembling  a 

reliable  analysis  system  to  perform      auto-  and 
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cross-^brrela  From  the  start,  it 

was  realized  that  once  iere  constructed,  the 

problem  of  attain)  letector  efficiencies 
would  have  to  be 

Detector  efficiency  is  f  letector  size 

and  neutron  sensitivity „  of  all          pi  -<dlly  available, 

the  Westlnghouse  rj  iber  provided  the 

-ill 
best  combination  of  .  Ivity  amperes/ 

neutron/cm  /sec)  and  nsitive  area),  -As 

such,  this  dete  gh  (in  the 

hohlraum  and  lattice  .  measure  the  spectral 

density^  however  1  First,  it  provided 

a  noise  signal  compeu  use  which  would  be 

studied  in  the  lattice  lax         ius  provided  us  with  inf or<= 

matlon  on  the  design  s  analysis  system „ 

Second,,  any  very  large  i  the  hohlraum  which  might 

be  present  due  to  spatial  effects         11  be  detected* 

Thus,  both  auto-  and  •  ,..,  were  made 

in  the  lattice  to  insure  that  a  reliable  system  would  be 

available  in  the  future  0  Concurrently,  a  separate  study  was 

being  made  of  methods  for  in   ising  det*     efficiencies 

to  the  levels  needed  A      * 

Chapter  2  describes  the  he  Initial  phase 

;embl;v  Chapters  3 

I  th 

Chapter  5 
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CHAPTER  2 

ASSEMBLY  OF  EQUIPMENT  AND  PRELIMINARY 
HOHLRAUM  MEASUREMENTS 

2ol  Assembly  of  Equipment 

Several  methods  of  analyzing  noise  signals  have  been 
developed,  often  classified  as  to  whether  the  analysis  is 
done  concurrently  with  reactor  operation,  or  done  afterward 

by  analyzing  either  an  electronically  or  mechanically  (pen- 

.  (q) 

recorder)  recorded  signal,1    Whether  the  analysis  is  done 

digitally  or  continuously,  further  specifies  the  method  of 
analysis  to  be  used.  The  choice  of  method  of  course  depends 
upon  financial  and/or  time  limitations,  and  also  upon  the 
degree  of  accuracy  desired . 

For  present  purposes,  it  was  decided  that  the  use  of 
analog  computing  equipment  would  provide  the  most  flexible 
system  of  analysis*  Not  only  would  it  provide  a  real-time 
method  of  continuous  analysis  but  it  would  also  be  adaptable 
to  other  uses  such  as  reactor  dynamics  simulation.  For  this 
reason  one  George  A0  Philbrick  MU/DV  Duplex  Multiplier/Divider 
and  three  G.A0P0  Model  HK  Operational  Manifolds  (each  con- 
taining five  chopper-stabilized  operational  amplifiers),  were 
acquired.  The  latter  afford  the  capability  of  performing  a 
total  of  fifteen  separate  linear  operations  on  one  or  more 
signals,  whereas  the  former  performs  the  non-linear  operations 
of  multiplying  or  dividing,   (See  Appendix  D  for  a  brief 
discussion  of  the  characteristics  of  operational  amplifiers 
and  a  description  of  the  basic  oi       \ised  in  performing 
the  linear  operations „ )  The  method  oi       is  used  was  to 
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take  the  Fourier  integral  .  function 

with  zero  delay  time.   A  )  the  ci]    ry   is 

shown  in  Figure  2.10 

After  using  the  he  technique  of  multiplying 

the  detector  noise  by  a  sine  wave  of  known  frequency,  a  low 
frequency  band-pass  filter  of  fixed  bandwith     is  used  to 
pass  only  those  noise  components  whose  frequency  differs  from 
that  of  the  oscillator  by  less  than  the  filter  bandwidth. 
The  resulting  si^na,  on  by 

the  integration  time  gives  fcl  spectral 

density  at  a  frequency  cu>  -  Aw,  In  a  bandwidth  A^  0 
Scanning  the  spectrum  Is  accomplished  by  varying  the  oscilla- 
tor frequency  to      over  the  de 

2 . 2  Pre 1 iminary  Experiment at I on 

Initial  assembly  of  the  equipment  and  e    mentation 
were  done  by  the  author  togeth* 

As  stated  earlier^  this  first  phas  was  I        wo  purposes « 
First  to  determine  if  the  analyzing  was  performing 

correctly,  and  secondly  to  de  thermal  column 

and  hohlraum  served  as  a  true  white  noise  source  of  neutrons. 
This  latter  question  is  impoi  seen  by 

Eq.  (l015),  determination  of  the  lattice  transfer  function 
requires  foreknowledge  of  the  spe      density  ©f  the  source. 
Only  when  the  external  source  is  itself  random  may  the 
simplified  pile  noise  formulation  be  used, * 

For  these  measurements,  a  We       se  Model  6377  Ioniza- 
tion chamber,  was  placed  in  the  h  irough  a  hole  bored 
through  the  graphite  wall  inside  port 
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15 
Preliminary  measurements  showed  that  the  total  cLe0 

detector  output  was  175  ixamps  with  a  noise  component  of 
about  1  uamp„   Since  it  was  desirable  to  have  a  voltage 
input  to  the  vacuum-tube  equipment  on  the  order  of  one-half 
volt  or  larger,  the  need  for  a  high  gain  preamplifier  became 
apparent c  To  avoid  excessive  60  cycle  pickup  and  to  intro- 
duce as  little  instrument  noise  as  possible,  it  was  decided 
to  use  an  Analog  Device  Model  106  Solid  State  Aoco  opera- 
tional amplifier.,  powered  by  four  F4BP  6  volt  batteries  „  In 
addition  to  meeting  technical  requirement ss   the  amplifier 
has  the  added  advantage  of  economy  costing  $26 0  For  similar 
reasons,  the  detector  high  voltage  was  supplied  by  10  67„5  V 
batteries .  To  avoid  saturating  the  preamplifier  with  the 
do  Co  level  of  the  incoming  signal  while  still  obtaining  maxi- 
mum amplification  of  the  noise  component ,  a  buckling  voltage 
controlled  by  a  potentiometer  was  inserted  just  before  the 
input o      Also,  since  there  was  a  large  low  frequency  variation 
in  the  neutron  flux,,  apparently  due  to  control  rod  motion  in 
the  reactor*  a  high  pass  network  (Rg  and  C2)  with  a  break- 
point at  205  cps  was  introduced  (see  Figure  202)0  Thus  the 
preamplifier  response  iss 

~2R2R1^  +  2^2 C2S^         *  (o  ,\ 

e     -  — — —     i    —  i  i^oi; 

lip     L»pS 

^2R2   +  Rl^1   +  2L+T' 
which  for  large  s  =  J«o    becomes 

•r~-  =  -  R»   o  I  -  °  2 ) 

Although  a  high  frequency  roll-off  due  to  finite  space  charge 
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collection  times  was  observed,  the  preamp-detector  frequency- 
response  was  constant  up  to  5*000  cps  for  a  detector  voltage 
of  720  V, 

The  noise  amplitude  at  the  preamp  output  was  about  one- 
half  volt,  which  was  sufficiently  large  compared  to  the 
0o01  V  self -noise  of  the  preamp .  Two  stage  amplification  with 
a  total  gain  of  200  provided  a  100  V  peak-to-peak  signal  at 
the  input  to  the  heterodyning  multiplier.  Here,  the  signal 
is  multiplied  by  a  100  V  peak-to-peak  sine  wave  of  known 
frequency.,   Since  the  output  of  the  multiplier  is  scaled  by 
100  V,  the  resultant  signal  is  100  V  p-p0  The  combination  of 
the  two  filters  produced  a  0088  cps  pass-band  centered  at 
0,128  cps.  The  inclusion  of  an  amplifier  with  a  gain  of  20 
between  filters  maintained  the  signal  at  a  sufficiently  high 
level  to  make  the  effects  of  component  drift  at  this  stage 
negligible o   After  squaring,  the  filtered  signal  was  integra- 
ted using  an  averaging  or  "leaky"  integrator 5  with  a  time 
constant  of  100  seconds „  The  circuit  at  this  point  is  shown 
in  Figure  2,3°  Results  of  the  preliminary  measurements  in 
the  hohlraum  are  presented  in  Figure  204o  Two  features  are 
important o  The  first,  already  described,  is  the  high 
frequency  roll-off  as  a  function  of  detector  voltage.  The 
second  is  the  low  frequency  erratic  peaking  in  amplitude . 

To  investigate  the  latter,  the  effect  of  the  reactor 
source  neutrons  was  considered.   Since  these  neutrons  orig- 
inate in  the  critical  MIT  Reactor,  a  low  frequency  breakpoint 
in  the  reactor  neutron  spectral  dens"     -;ld  be  expected  at 
^  =  p/_/  .  This  value  is  approximately  0 00757*0015  for  the 
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19 
MIT  Reactor*  resulting  In  a  break  at  f  »*-/2ir  -  0o796  cps„ 

In  order  to  determine  the  amount  of  correlation  intro- 
duced into  the  lattice  by  these  source  neutrons^  one  can 
assume  that  the  lattice  is  a  detector  and  apply  Eq0  (l„l6)0 
To  do  this^  the  efficiency  of  the  lattice  for  reactor  neu- 
trons must  be  estimatedo   Less  than  0o02#  of  the  reactor 
neutrons  pass  the  hohlraum0 *  '  This  value  may  be  taken  as 
an  upper  limit  for  the  efficiency  of  the  lattice.  Using 
Eq.  (lcl.6)  for  the  case  of  a  critical  reactor  (kfc£  l)9   with 
a  =  4o8  sec  9   it  is  seen  that  above  about  5  cps^  correlation 
of  neutrons  in  the  lattice  should  be  negligible „   Since  the 
frequency  range  of  interest  in  the  lattice  is  about  3^000 
cycles 5  the  effect  of  the  low  frequency  correlation  may  be 
ignored o  Thus  the  hohlraum  may  be  considered  to  be  a  white 

noise  source  of  neutrons  for  the  lattice c 

-tVie 
For  this  reason,,  tb4«  behavior  of  the  spectral  density 

curve  of  Figure  2„4  below  100  cps  caused  some  concern^ 

2 °  3  First  Phase  Conclusions 

In  order  to  check  on  the  reliability  of  the  system., 
measurements  of  the  spectral  density  were  taken  several  times 
on  different  days.  Above  50  cps*  the  results  were  found  to 
be  reproduceable  within  the  predicted  experimental  error  (see 
Appendix  D)„   Scattering  of  data  points  below  50  cps  cast 
doubt  upon  the  reliability  of  the  system  at  low  frequencies 0 

Preliminary  investigations  were  concluded  at  this  point 
with  the  next  objectives  being  tog 

a)  attempt  to  improve  detector  efr  Disney 
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b)  investigate  more  thoroughly  the  low  frequency 
.  response  of  the  system 

c)  investigate  frequency  response  in  the  lattice 
itself. 

The  problems  of  developing  a  detector  with  increased  effici- 
ency for  lattice  noise  analysis  is  the  subject  of  another  SM 
thesis1   '  by  McG.  Johnson*  currently  underway .   The  next 
chapter  describes  work  done  by  the  present  author  on  the 
items  (b)  and  (c)„ 
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CHAPTER  3 

LOW  FREQUENCY  TESTS  AND  AUTOCORRELATION  MEASUREMENTS 

Section  2.3  described  the  erratic  behavior  of  the  analyzer 
at  low  frequencies o   It  was  necessary  to  determine  the  cause  of 
this  behavior  before  proceeding  to  lattice  autocorrelation 
measurements,, 

3»1  Introduction 

Up  to  this  pointy  the  oscillator  signal  had  been  provided 
by  first  a  General  Radio  Unit  R~C  Oscillator  Type'1210-C  and 
later  a  Hewlett  Packard  Test  Oscillator ,  Model  65OA*  the  latter 
providing  the  signals  of  lowest  frequency  down  to  10  cps„   In 
order  to  look  more  closely  at  the  low  frequency  response^,  a 
Waveforms  Model  403B  oscillator  was  used,  providing  signals  as 
low  as  1  cpSo  The  filters  were  also  modified  as  shown  in 
Figure  3»L 

This  time  the  spectral  density  of  the  hohlraum  noise  was 
investigated  in  the  range  1  to  200  cps5  with  the  results  shown 
in  Figure  3<,20   As  can  be  seen*  results  were  again  erratic,, 
with  a  general  trend  of  higher  amplitudes  at  lower  frequencies. 

The  erratic  low  frequency  response  suggested  the  possibil- 
ity of  component  drift *  while  the  increase  in  amplitude  at 
low  frequencies  suggested  that  the  d0c0  portion  of  the  detec- 
tor signal  was  not  being  completely  cancelled*  allowing  the 
carrier  wave  to  be  passed  to  the  filters.  Referring  to  the 
latter  problem*  note  that  a  one-half  volt  d0c„  signal  at  the 
noise  input  to  the  multiplier^  when  multiplied  by  a  100  V 
peak-to-peak  carrier  signal^  of  f  reque>       will  result  in 
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25 
a  one-half  volt  peak-to-peak  sine  wave  at  the  multiplier 
output o   At  high  frequencies  (such  that  uj      is  much  greater 
than  the  center  frequency  of  the  filter)  this  component  is 
filtered  out  leaving  only  the  signal  due  to  the  a»c0  (noise) 
portion  of  the  signal „  However ,  since  the  filters  are  not 
square  (effectively  24  db /octave  attenuation  as  seen  at  the 
averaging  circuit  input),  at  low  frequencies  this  sine  wave 
would  not  be  attenuated  enough  to  be  negligible.  Total  cancel- 
lation of  the  do  Co  level  of  the  signal  was  insured  by  adding 
a  blocking  capacitor  in  series  with  the  first  amplifier., 

Multiplier  drift  (rated  at  +1$)  which  was  measured  to 
be  as  much  as  2  V  in  one  hour,  accounted  for  the  erratic 
behavior  at  low  frequencies „     This  drift  was  minimized  by 
adding  a  ventilating  fan  to  provide  better  temperature  stabil- 
ity? however,  drifts  of  up  to  0=75  V  per  hour  were  still 
encountered.,  Figure  3o3  shows  how  these  measures  served  to 
reduce  these  problems „ 

3o2  Simulation  of  Lattice  Measurement 

To  Insure  that  the  equipment  was  indeed  operating 
properly,  it  was  decided  to  simulate  a  lattice  measurement 
by  adding  a  capacitor  to  the  preamplifier,  in  parallel  with 
the  feedback  resistora  Neglecting  the  low  frequency  filter 
network,  the  preamplifier  has  a  gains 

JL  '  1  +  R-C-.B   °  {3o±} 

n       ir 

With  R^  -  S10K  (1#)  and  Cf  =  ured)  this  results 

in  a  simulated  breakpoint  at? 
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2iR^C~  =  314  ±  4  cps 


At  the  same  time,  the  filters  were  modified  to  widen  the 
pass  band,  which  would  allow  shorter  integration  times  (see 
Appendix  D).  The  filter  modifications  are  shown  in  Fig„  3„40 

Results  of  the  simulated  lattice  measurement  are  shown 
in  Fig0  3o5o  The  breakpoint  occurs  at  that  frequency  corres- 
ponding to  a  6  db  decrease  in  the  squared  signal  and  is 
estimated  to  be  320  +  20  cps5  which  is  in  good  agreement  with 
the  predicted  vali. 

3 o 3  Lattice  Autocorrelation  Me a sur ement s 

Having  been  satisfied  that  the  system  was  now  functioning 

properly ,,  attention  was  turned  to  making  a  test  run  in  a 

D00-moderated,  uranium -fueled  lattice.   Estimates  of  k'   and 
<=  eff 

a  for  the  lattice  installed  in  the  exponential  facility  during 
this  thesis  work  are  developed  in  Appendix  C,  indicating  that 
a  detector  efficiency  of  around  30$  would  be  required  to 
enable  determination  of  a  for  this  particular  lattice.,  None- 
theless, since  it  would  be  several  months  at  least  until  a 
more  reactive  lattice  was  scheduled  for  installation,  it  was 
decided  to  make  a  test  run  to  insure  that  the  equipment  would 
be  ready  for  use0  These  measurements  were  made  by  suspending 
the  6377  boronlined  detector  in  a  lead-weighted  aluminum  tube 
in  the  center  of  the  lattice  core0   In  order  to  do  this,  the 
seven  center  fuel  rods  had  to  be  removed  from  the  211  rod 
lattice 0  Since  the  multiplication  factor  of  the  assembly  w 
already  so  low  as  to  preclude  measm-     of  kinetic  effects, 
the  small  additional  decrease  in  react iv;?t;"  caused  by  removing 
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the  rods  was  not  considered  significant ,  First  attempts  at 

taking  measurements  revealed  that  the  75  foot  cable  from  the 
detector  to  the  preamp  caused  excessive  oscillation  in  the 
preamplifier  <>   Although  this  oscillation  could  be  eliminated 
by  use  of  an  RC  compensating  circuit,,  the  signal  strength  was 
then  so  small  as  to  be  comparable  to  the  preamp  self -noise . 
Therefore,,  the  preamplifier  was  placed  on  the  lattice  top* 
eliminating  the  need  for  both  the  long  input  cables  and  the 
compensating  networks   Such  an  arrangement  produced  lattice 
noise  signals  five  times  greater  than  the  preamp  self  noise „ 
With  the  equipment  set  up  in  this  manner 5  the  results  shown 
in  Figure  3 .6  were  obtained „   Again,  the  white-noise  spectrum 
due  too  the  low  detector  efficiency  is  evident „ 

No  further  work  was  performed  on  autocorrelation  measure- 
ments because  the  cross-correlation  technique  to  be  described 
In  the  next  chapter  appeared  to  hold  more  promise 0 
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CHAPTER  4 

CROSS-CORRELATION  MEASUREMENTS 

4ol     Theory 

As  mentioned  in  Chapter  1^  cross-correlating  the 
signals  should  decrease  the  detector  efficiencies  needed 
for  noise  measurements .  The  cross-correlation  function 
is  defined  as  the  expectation  value  of  two  independent 
signals s 

Although  0-,   and  #L  consist  of  noise  signals  super- 
imposed on  doCo  levels,  we  need  consider  only  the  noises 

$fx{t)   =u1(t)  +  c1(t)  (4,2) 

02(t)   =U2(t)  +  C2(t)  ,  C403) 

where  U  and  C  are  uncorrelated  and  correlated  noise , 
respectively o 

The  cross-correlation  function  of  these  two  signals 
is  given  bys 

E^CtJUgCt)]    +  EJU^tJC^t)]   +  EJUgftJC^t)]  +  EJO^tjGgft;) 


Since  the  expectation  value  of  the  multiple  of  two 
uncorrelated  functions  is  zero,  all  uncorrelated  terms 
disappear  ands 

B^(*V2<*J]  =  EJc^tJC^t)]  .  (M) 

Thus,  corresponding  to  EqQ  l„l6  for  auto-correlation 
experiments,  one  can  write  for  cross-correlation  experiments  § 
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<lIc(p)|>2=2€16sP0Q2 


where  now  €  ±   and  £ 2  refer  to  the  efficiencies  of  the 
separate  detectors .  This  would  seem  to  imply  that  measure- 
ments could  be  made  independent  of  detector  efficiency,  as 
long  as  the  analyzer  output  signal  is  sufficiently  large c 
However,,  as  stated  earlier,  Sei^fritz,  et  ale  found  that 
cross-correlation  decreased  their  efficiency  needs  only  by  a 
factor  of  20 „  This  is  due  to  the  fact  that  the  cross-terms 
in  Eq0  (4»4)  are  not  totally  uncorrected <>  A   factor  of  20 
reduction  in  efficiency  requirements ,  however^  is  a  highly 
desirable,  if  not  essential,  improvement  if  one  is  to  make 
noise  measurements  in  subcritical  systems „ 

ho2     Equipment 

The  cross-correlation  measurements  required  the  addit- 
ional equipment  shown  within  the  dotted,  lines  in  Figo  4ol, 
consisting  of  one  additional  multiplier  (G,A,P.  +  25  V  MtJ/DV 
Duplex  Multiplier/Divider  )<>  6  operational  amplifiers, 
another  preamplifier,  and  an  additional  Westinghouse  Type 
6377  boron-lined  ionization  chamber  * 

The  two  detectors  were  suspended  in  the  high-flux 
region  of  the  lattice  described  earlier,  side-by-side,  in 
separate  lead-weighted  aluminum  tubes „  Again  the  preampli- 
fiers were  placed  on  the  lattice  top,  as  close  to  the  detec- 
tors as  possible  to  avoid  preamp  oscillations «  Both  detec- 
tors were  operated  at  670  V° 
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4o3  Result's' 

Recall  that  for  autocorrelation  measurements  the  filtered 
signal  is  squared  and  then  averaged „  The  squaring  process 


naturally  results  in  a  large  non-geJsper  signal  *  whose  average 


'It     1  ^  X, 

\/Jf^            2 

v|Ic'  r     ewv-i) 

difll^       (v)2 

J  I  £  p 

will  also  then  be  non°3ej*©0  The  cross-correlation  measurement^ 
in  that  it  eliminates  the  uncorrelated  noise  signals*  should 
produce  a  signal  whose  magnitude  relative  to  the  autocorrela- 
tion signal  is  given  by  the  ratio  of  Eq0  (4»6)  to  Eq»  (1.16). 
For  very  low  efficiencies^,  this  ratio  is  merely  the  right-hand 
term  of  Eq<>  (lol6)  and  iss 

(4.6) 

For    £c?^10~'*   and  for  the  lattice  under  study  % 
<fJXQ<2>  .6 

</iar> 

Thus  since  the  signal  level  in  the  auto-correlation  measure- 

ments  was  dX   I  >-  100  V*  the  cross-correlation  output  signal 

a 

must  be  <IIQI2>  -  3  x  10"  V.  This  signal  is  so  low  as  to  be 
overridden  by  equipment  drift  „  While  observing  this  drop  in 
net  output  confirms  predictions  as  to  the  effect  of  converting 
from  auto-  to  cross-correlation*  the  low  signal  made  it 
impossible  to  extract  a  spectrum „  Even  with  cross-correlation 
techniques*  therefore*  higher  lattice  multiplications  and 
detector  efficiencies  will  be  required „ 

Experimentation  was  concluded  at  this  point ,  Conclu- 
sions and  recommendations  follow  in     ^r  5° 


36 

CHAPTER 

SUMMARY,  CONCLUSIONS ,  AND  RECOMMENDATIONS 


A  reliable  system  for  making  auto-  and  cross~power  noise 
spectral  density  measurements  on  a  subcritical  nuclear 
assembly  was  constructed „  The  system  was  demonstrated 
capable  of  providing  signal  analysis  over  a  frequency  range 
of  1  to  5*000  cpso 

The  assembled  system  was  tested  in  several  applications? 

a)  a  signal  simulating  a  reactor  noise  spectrum 
was  analyzed  to  demonstrate  successful  system 
performances  under  optimum  operating  condi= 
tlons  (Chapter  3$   Section  2) 0  The  measured 
breakpoint  frequency  320  +20  ops  agreed  quite 
well  with  the  known  input  of  314  +  4  epso 

b)  The  MIT  exponential  facility  hohlraum  was 
studied  using  the  autocorrelation  technique  „ 
None  of  the  results  (Chapter  3<>  Section  1)  gave 
reason  to  doubt  the  assumption  that  the  hohlraum 
provides  a  white  source  of  neutrons  to  the  test 
lattices . 

c)  Both  auto-  and  cross~correlation  measurements 
were  made  in  a  lattice  „  In  both  cases  it  was 
shown  that  lattice  multiplication  and  detec° 
tion  efficiencies  must  be  increased  before 
noise  analysis  can  be  used  t ?  extract  neutron 
decay  constants  (Chapters  3  and  4) « 
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Since  detector  efficiency  requirements  vary  inversely 
as  the  snuare  of  the  multiplication  factor  of  the  system, 
future  studies  should  concentrate  on  lattices  with  k  rr   as 
close  to  the  present  lattice  safeguard  limit  of  k  f  f  =  0„9 
as  possible ,  Since  cross -correlation  of  two  signals 
significantly  reduces  the  efficiency  requirements,  this 
method  provides  the  greatest  promise  for  successful  analysis 
at  reasonable  efficiency  levels „ 

3°2  he  c  omme  nda  t  ion  s 

In  order  to  attain  the  required  detector  efficiencies, 
the  utility  of  gamma  detectors  in  the  lattice  annulus 
should  be  fully  explored-  Such  detectors  would  not  disturb 
the  multiplying  system,  but  would  detect  both  prompt  fission 
product  gammas,  which  escaped  from  the  lattice  core,  and 
those  gammas  produced  by  neutron  capture  in  the  cadmium 
sheath  surrounding  the  lattice  „  Alternatively  it  may  be 
possible  to  develop  a  detection  system  which  detects  the 
Cerenkov  radiation  produced  in  the  moderating  DgO.  Further, 
the  2$>   enriched  UOp  lattices,  scheduled  for  study  beginning 
in  February  1967,  should  provide  the  capability  of  attain- 
ing k  ff  ^0o9,  which  would  also  reduce  detector  efficiency 
requirements  considerably « 

Finally,  in  future  work,  where  large  volume  detectors 
must  be  used,  the  space  dependence  of  reactor  noise  measure- 
ments should  be  carefully  considered. 
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APPENDIX  A 
NOISE-EQUIVALENT  SOURCE  FOR  A  SUBCRITICAL  ASSEMBLY 

Neutron  processes  in  a  multiplying  system  are  analogous 
to  the  random  flow  of  electrons  in  a  diode „  Using  this 
similarity ,  Cohn  has  adapted  the  diode  noise  equation  of 
Schottky  to  the  noise  in  a  critical  reactor 0C10)  By  simple 
modifications  it  is  possible  to  extend  this  treatment  to  the 
case  of  a  subcritical  system <, 

Cohn's  adaptation  of  the  Schottky  formula  iss 

<|S  !2>  =22  q  2i  ,  (Ad) 

i 

2 

where s  CIS   I    >  is  the  noise ^equivalent  source; 

q.  -  net  number  of  neutrons  produced  in  the  nuclear 

process  of  type  i;  and 
m.  =  average  number  of  processes  of  type  i  occurring 

per  second  in  the  system . 
Table  A.l  lists  the  processes^,  their  rate  of  occurrence, 
and  the  net  number  of  neutrons  produced  in  that  process „ 

TABLE  Aol 

CONTRIBUTIONS  TO  PILE  NOISE  SOURCE 

Average  Rate        Net  Number  of 
Nature  of  Process      of  Occurrence      Neutrons  Produced 


n   A 


Nonproductive  leak- 
age or  absorption       <£       + 

Fission  giving  rise      n   F   p  N=1 

to  N  prompt  neutrons 

External  source  provide 
ing  neutrons  to  main- 
tain the  pile 


n 

F 

PN 

JiA 

+  F 

7(: 

I  -  kp) 
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In  this  formulation  delayed  neutron  effects  have  been  negle< 
tedo  Here 

A  =  macroscopic  cross-section  for  all  non-productive 

processes  including  leakage , 
F  =  macroscopic  cross-section  for  fission* 
n  =  total  number  of  neutrons  in  the  system* 
Jl  =  prompt  neutron  lifetime  in  sec,  and 
A  =J?/k  =  neutron  generation  time  in  sec  a 
A  and  F  are  subject  to  the  neutron  balance  condition 


VF 


\M  ad) 


A  4-  F   ""p 

where  &   is  the  average  number  of  neutrons  produced  per 
fission o  P  is  the  probability  that  N  prompt  neutrons  will 
be  produced  in  any  one  fission*  subject  to  the  conditions 

oo 

£  PM  =  1  (Ao3) 


CD 

S  N  PM  =  (1= 
N=l   w 


(A  .4) 


N=l   W 
Substitution  into  Equation  (Aid)  yields 


(A*5) 


<,so|2>=^T^ 


A+kF  S  (N~l)2PM  +  (l-k)   o  (Ao6) 


Applying  conditions  (Ao3)*  (Ao4)  and  (A„5)  yields 


<>S0<2>  -f- ; 


•**       x. 


4-(?-2P+l] 


(A  ,7) 


The  first  term  in  square  brackets  is  t    ->rm  due  to 
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nonproductive  (e0g0,  (11,7)  capture  and  leakage)  processes, 
the  second  is  the  term  due  to  fission  processes,  the  third 
is  the  term  due  to  the  external  source „  Note  that  the  exter- 
nal  source  term  has  been  assumed  to  be  frequency  independent, 
and  that  the  entire  noise  equivalent  source  is  thus  a  white 
noise  source 0  Since  the  fission  rate  iss 
nk_ 

*o  --=5  »  (A°8) 


EqQ  (Ao7)  becomes 


5> 


k 


2 


k^  +  -£— (v 
P    "v 


2V   +  1 


(A.9) 


The  source  neutrons  are  then  correlated  by  the  multiplying 
process •  The  correlated  component  of  the  noise  spectrum  is 
obtained  by  multiplying  equation  (A»9)  by  the  nuclear  system 
transfer  functions 

i2^    ■„./*.%  i2  •,„  ,2 


</n(f)r>  =  iT(f)r  <is0r> 


2)JF 


o 


k 


P 


2. 


k        k' 

-3-  -  K  +  -Av*-25  +y  ) 

v        p   v 


OJ 


1 

2  \JL 


1-k. 


(A  .10) 


where  a  =    P  is  the  prompt  neutron  decay  constant  0 

The  average  current  I  in  a  detector  due  to  all  fissions 
detected  will  beg 


nk 
I  -  €q  -=4  =£0Fo 


(Aoll) 


where  Q  is  the  charge  collected  per  fission  detected „     Super- 
imposed  on  this  d0c0   level  are  noise  components  due  to  the 
correlated  neutrons  and  the  randoranec        f  the  collection 
process .     The  former  is  given  by  Eqs „   (A .10)  and  (Aoll) 


p  (>)3T 
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(A  .12) 


262Q2k 


£  p 


{v)f 


2  -  J£  -  k     +  i*-^  -  2V+  l)U^r  • 
v         P        y  ~L>d+ad 


(a  as) 

The  latter  can  be  obtained  by  using  the  Schottky  formulas 


<|IC|2>  =  2Q2€FQ 


(A  .14) 

For  £<  <  1,  these  two  terras  are  first  order  uncorrelated 
and  are  added  in  quadrature  to  give 

€k  -    ..  ,.2 


</ia(^)r>  =  2q^p0 


ck    r       k  k 

1  +  Jib  -  JB.  «  k     +  ^-(^=27  +1) 

I  STL         »  P         " 


(A.15) 
This  result  enables  estimation  of  detector  efficiencies 
needed  to  extract  information  from  noise  signals  as  shown 
in  Chapter  1,  Section  4. 

Note  that  throughout ,  it  has  been  assumed  that  the 
analysis  equipment  did  not  introduce  any  spurious  frequency 
dependence .  Such  an  assumption  is  justified  by  the  "white" 
noise  spectrum  displayed  in  Figure  3o3o 


APPENDIX  B 
OPERATIONAL  AMPLIFIERS  AND  ANALOG  CIRCUITS 
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An  operational  amplifier  is  a  high  gain  amplifier  with 
certain  input  and  output  impedances  as  shown  in  Figure  bu1o 
Analysis  of  the  circuit  shows  thatg 


ei  *T 


1+I<Z7+1> 


(Boi) 


where  Zf  =  feedback  impedance 

Z-  =  input  Impedance 

A  =  open  loop  gain* 

Zf 
Note  that  if  A  >>(—-+  1),  this  reduces  to§ 


(B.2) 


ei = "  zi  • 

Since  Zf  and  Z.  are  complex  quantities,  various  combinations 
of  resistances  and  capacitances  may  be  used  to  perform  mathe- 
matical operations o  Table  Bol  Illustrates  the  circuits 
which  are  used  to  perform  the  linear  operations  in  the 
analyzing  equipment „ 


Figo  Bol  Operational  Amplifier  Crr-cuit 
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APPENDIX  C 

ESTIMATION  OP  keff  AND  a 

Estimates  of  k   and  a  can  be  made  by  using  the 
THERM0S  computer  code 9^      '     This  integral-transport  theory 
code,  based  on  the  Wigner-Seitz  unit  cell  models  permits 
the  calculation  of  § 

f,   the  thermal  utilization 

*7   the  average  number  of  neutrons  produced 
per  thermal  neutron  absorbed  in  the  fuel 

2   the  average  macroscopic  absorption  cross- 
section  in  the  cell 


of 


Combining  these  values  with  reasonable  estimates 

t  -  age  to  thermal 

£  =  fast  fission  factor 

p  =  resonance  escape  probability 

D  =  diffusion  constant 

B  =  geometric  buckling, 
g 

k   can  be  estimated  bys 

(VKpf)exp(-B  2t) 
k  \     r?  ,  (Col) 


00      1  +  L2^2 

6 


2 


where  L2  =  —  =  (diffusion  length) 

The  lattice  investigated  in  the  present  studies  had 
the  following  specifications s 
Pitch  2,25  in. 

Enrichment  0a$W   U235 

Fuel  Density       18*90  gm/cc 
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D20  Purity  99  .4  MPC 

The  unit  cell  dimensions  (outer  radius  minus  inner  radius) 

ares 

Fuel  0.49149  cm 

Air  gap  0,00762  cm 

Cladding  0.08636  cm 

Moderator  2 .41491  cm 

Using  this  data,  and  the  prooe&ures  outlined  in  reference 

(14),  the  following  estimates  were  derived  from  the  THERMOS 

outputs 

2  =  0.007583  cm"1 
f  =0.969 
y(     =  1.447 

For  this  lattice,  good  estimates  for  the  other  quanti- 
ties^ '  ares 

D  -  8»75  cm 
p  *  0.92 

€  -  1 oOl 

I  v  115.4  cm2 

B„  =  0.0031  cm"2 
g 


These 

values 

result  ins 

koo 

-j 

1.30 

keff 

*-\j 

0.65 

1. 

JL 

*w 

4.35  x 

10"     sec 

a 

1  ~ke: 

^£.^805  sec4 

Consulting  Figure   1.2  one   sees 

in  this  lattice, 
required  efficiencies  for  successful  noi&e  analysis  experi- 
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merits  would  bes 

for  autocorrelations       30# 
for  cross-correlation     1.5#  „ 
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APPENDIX  D 

COMMENTS  ON  OPERATION  OF  EQUIPMENT 

Dol  Equipment 

The  equipment  assembled  to  perform  the  auto  spectral 
analysis  consists  of  the  following  (see  Figs0  D0i  and  D02)8 

1  George  Ac  Philbrick  (GAP)  MU/DV  Duplex  Multiplier 
Divider  -  +  100  V  input -out put ,  Ser*  0297 

2  GAP  Model  HK  Operational  Manifolds*  Ser0  O69  and  068 
8  GAP  K2-XA  Operational  Amplifiers 

8  GAP  K2-P  Stabilizing  Amplifiers 

8  GAP  2-Unit  K-l  Modular  Assembly  Units 

1  Battery  Powered  A »C „  Pre -amplifier  utilizing  Analog 
Devices  Model  106  Operational  Amplifier 

2  Associated  Specialties  Model  2  Power  Supply*  Ser* 
748  and  749 

1  Westinghouse  Model  6377  Ionization  Chamber 

10  NOo  467  67.5  Volt  Batteries 

1  Waveforms  Model  403B  Sine-Square  Wave  Generator, 
Ser0  7997 

1  Esterline  Angus  Graphic  D0C0  Ammeter  Model  AW* 
Ser0  85339 

Each  linear  operation  is  performed  using  one  K2-Xa  and 
one  K2-P*  which  are  plugged  side-by-side  into  the  top  of  the 
manifold o  Each  manifold  can  accommodate  5  pairs  of  these 
amplifiers .  The  K-l  Modular  Assemblies  are  labeled  to  design 
nate  which  operation  is  performed „  These  assemblies  are 
plugged  Into  the  face  of  the  HK  manifold  so  as  to  complete 
the  circuitry  for  each  chopper-stabilised  operational  ampli- 
fier o     Each  K-l  Assembly  is  labeled  as  to  input,  output,  and 
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Preamplifier 

Detector 
Power  Supply 


Multiplier 


K  Manifold 
K2-P  and  K2-XA   Pairs 
and   Kl  Modules 


— Waveforms  403B 
Oscillator 

Esterline  Angus 
-^Recorder 


Pig.  D.l  Auto-Spectral  Density  Analysis  Equipment 
(Connecting  Cables  not  Shown) 
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Preamplifier 

Detector 
Power  Supply 


Multiplier 


K  Manifold 
K2-P  and  K2-XA  Pairs 
and  Kl  Modules 


— Waveforms  403B 
Oscillator 

.Esterline  Angus 
^■"Rec  order 


Fig.  D. 1  Auto-Spectral  Density  Analysis  Equipment 
(Connecting  Cables  not  Shown) 
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Figo  D02  Photograph  Showing  Filter.  Module  (Left 
Preamplifier  (Right ) 
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ground  connections „  The  eight  operations  performed  are*  in 

order  (see  Fig0  4„l)s 

Ao  Gain  of  20  amplification  with  D„C.  blocked 

Bo  Gain  of  10  or  20  amplification  for  noise  signal 

Co  Variable  gain  (up  to  20)  amplification  for  oscillator 

signal 

Do  Band  pass  filter  which  may  be  switched  to  provide  either 

narrow  band(^003  cps)  or  wide  band  (^10  cps)  operation 

Eo  Gain  of  XXj  X5$   X10,  X?Q*  or  X100  switchable  amplifier 

F.  Second  filter  similar  to  D 

Go  Averaging  circuit 

The  time  constant  of  the  averaging  circuit  may  be  altered 
by  using  banana  plugs  to  insert  various  combinations  of 
parallel  RC  components  into  the  jacks  marked  FEEDBACK «» 

In  order  to  perform  the  cross  power  spectral  measurements, 
the  following  equipment  was  duplicated 

detector 

preamplifier 

detector  battery  power  supply 

operational  amplifiers  and  modular  units  to  perform 
operations  A*  B9   D,  E^  and  F 

A  +  25  V  GAP  MU/DV  Multiplier  was  borrowed  from  the  MIT 
Mechanical  Engineering  Department  to  serve  as  the  correlating 
multiplier o  A  GAP  Model  R-300  Power  Supply  was  also  borrowed 
to  meet  the  increased  power  requirements „ 

D<>2  Equipment  Operation 

In  order  to  insure  proper  opera g     the  following  had 
to  be  observed g 


51 

a)  Keep  equipment  well  ventilated   (e0g°  by  fan) 

b)  Maintain  largest  signal  possible  throughout  with- 

out exceeding  rated  signal  input  limitations 

c)  Avoid  ground  loops 

d)  Minimize  d.c0  amplification  at  preamp  by  varying 

resistances  in  blocking  circuit 

e)  Avoid  preamp  oscillations  by  keeping  input  cables 

to  a  minimum  length 
These  comments  supplement  the  individual  manufacturer °s 
manuals ,  which  have  been  included  in  a  booklet  Noise  Spectrum 
Analyzer^  compiled  by  the  author e 
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APPENDIX  E 

PROBABLE  ERROR  ESTIMATION 

Eol  Predicted  Error  In  Amplitude 

The  spectral  measurements  are  made  using  an  averaging  cir- 
cuit with  a  time  constant  X  =  RfCf  where  Rf  and  Cf  are  the 
feedback  resistance  and  capacitance  in  the  averaging  cir- 
cuit .  Using  formulas  developed  by  van  der  Ziel'  ^'  and 

f  -\f.\ 
Bendat*   ',  the  relative  standard  deviation  in  cross~power 

spectral  density  measurements  is  8 

"2  +  Nx/S  +  Ng/S  +  (N^SMNg/s) 


*S 


(Eol) 


irAf£ 
where  S  =  strength  of  signal  of  interest , 

N-,  p  =  strength  of  unde sired  signals  in  channels  1  and  29 
Af  =  filter  bandwidth  ,„  in  cps#  and 

£"  m  averaging  circuit  time  constant  0 
If  detector  efficiencies  are  great  enough  to  detect  a 
level  of  correlated  noise ,  S  and  N  would  refer  to  the  corre- 
lated and  uncorrelated  noise ,  respectively 0 

For  most  measurements  in  this  work,  detector  efficiencies 
were  so  low  as  to  permit  measurement  of  only  the  uncorrelated 
signal  level o  For  this  cases 


S     2 
S   Mr 


I1/2 

(E.2) 

Preliminary  hohlraum  measurements  were  made  with 

Af  =  0o88  sec"1 
2;  =  100     sec 
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Final  hohlraum  and  lattice  autocorrelation  measure- 

«1 
ments  were  made  with  Af  =10  sec 

T,  =  100  sec 

These  result  in  predicted  amplitude  errors  of  11  S$> 

and  2„5#j  respectively^,  and  are  considered  to  be  suffic« 

iently  small  for  this  work.  If  in  the  future  *  greater 

accuracy  is  considered  desirable ,  it  may  be  attained  by 

altering  the  feedback  impedance  in  the  averaging  circuit 

so  as  to  increase  "£  =  Rf^f  ° 

E.2  Fredicted  Error  in  Frequency 

Error  in  frequency  measurements  arise  from  the  finite 
bandwidth  of  the  filters  and  the  lack  of  accuracy  in  the 
oscillator  dial  readings  „  Frequency  resolution  due  to 
finite  bandwidth  is  +  Af ,  and  all  equipment  dial  readings 
are  rated  to  within  +2$., 

For  example j  using  a  filter  with  a  bandwidth  of  10  cps 
to  measure  the  spectral  density  near  a  frequency  of  300  cps 
would  result  in  an  uncertainty  in  frequency  of  +10  cps  due 
to  finite  bandwidth  and  +6  cps  due  to  dial  error .  This 
results  in  a  total  uncertainty  in  frequency  of  +12  cps,  or 
about  3=7#o  This  figure  could  be  reduced  by  decreasing 
the  filter  bandwidth „  As  can  be  seen,  however,  changes  in 
Af  affect  both  frequency  and  amplitude  uncertainty,  thus  a 
choice  of  Af  should  be  governed  by  desired  accuracy  in  both. 
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